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' Abstract 



Decoherence is studied in an attractive proposal for an actual implementation of a quantum computer 
, based on trapped ions. Emphasis is placed on the decoherence arising from the vibrational motion 

Q\ ' of the ions, which is compared with that due to spontaneous emission from excited states of the 

ions. The calculation is made tractable by exploiting the vast difference in time scales between 
the vibrational excitations and the intra-ionic electronic excitations. Since the latter are several 
orders of magnitude faster, an adiabatic approximation is used to integrate them out and find the 
£^ ' inclusive probability P(t) for the electronic state of the ions to evolve as it would in the absence 

of vibrational coupling, and the ions to evolve into any state whatsoever. The decoherence time is 
found at zero temperature and for any number of ions ./V in the computer. Comparison is made 
with the spontaneous emission decoherence, and the implications for how trap voltages and other 
parameters should be scaled with N are discussed. 

1 INTRODUCTION 

A quantum computer (QC) is a proposed machine that would make use of the superposition principle 
of quantum mechanics to do certain types of calculations with an unprcccndented degree of paral- 
lelism. The concept of a quantum computer goes back to the early 1980's.|i] || It is an outgrowth 
of the view that any computation is ultimately a physical process. This view has tended to get 
obscured in recent times, especially as our understanding of the universal mathematical properties 
of computation and associated issues in complexity and computational theory has grown. If one 
thinks about the computer as a machine, however, in which wheels must turn, or electrons must 
be transported between transisitors, and so on, then it is somewhat surprising, at least in hind- 
sight, that it has taken us so long to start thinking about the implications of the laws of quantum 
mechanics for the turning of the wheels, transport of charges etc. To avoid misunderstanding it 
should be emphasized that quantum mechanics already plays a vital role in the functioning of all 
modern electronic circuitry, in that it underlies the concepts of energy bands, the effects of doping 
on semiconductors, the theory of conductivity, the nature of electron and hole statistics, and many 
other aspects of semiconductor and device physics. The ultimate equations that we use to describe 
circuits, such as Ohm's law, or Shockley's equation, are deterministic and classical. Quantum me- 
chanics provides the microscopic basis for these equations, and fixes the parameters appearing in 
them. To use these equations, however, one only needs to know the meaning of the terms which 
appear in them, and the scope of their applicability. 
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In this introductory section we will describe what a quantum computer is. In Sec. 2, we shall 
describe the ion-trap quantum computer proposed by Cirac and Zoller||, and summarise the results 
of an earlier calculation jij of the decoherence time in this device. An approximate account of this 
calculation, less technical than that in Ref. [Q is given in an Appendix, along with some details about 
the structure of the ion array in the trap. Since many reviews of quantum computers now exist, Q 
this introductory section may well be redundant for many readers. Readers who are unfamiliar with 
the subject or with the terms used to describe the contents of the rest of the paper may, however, 
profit from it. 



1.1 What is a Quantum Computer? 

The key point is that the state of a classical computer is definite at all times. This is the prevailing 
rule in classical mechanics. In quantum mechanics by contrast, the state of a physical system is 
not definite, and a system does not have definite properties until they are measured. The well 
known uncertainty principle, which expresses the impossibility of simultaneously knowing both the 
momentum and the position of a particle with dcfinitcness, i.e., with zero uncertainty, is an expression 
of this fact. In general, therefore, the state of a quantum system must be expressed as a linear 
superposition of some basis states. Time evolution preserves this superposition property, and since 
it also conserves probability, must be described by a unitary operator, as long as the system is not 
externally disturbed. 

To understand how these ideas can be usefully applied to make a computer, it is best to take a 
simple example, that of evaluation of an integer valued function of an integer argument. The first 
point has to do with the correspondence between numbers and physical states. Thus a single bit 
(0 or 1) is represented by a physical system or register which can exist in two states, which are 
represented by |0) and (There are many examples of a two-state system of which we mention 
only two: the spin of a particle of spin-1/2, and the electronic states of an atom or ion. In the 
latter case, there are un general many more than two states, and the two-state description is only 
an approximation. One must be careful to avoid perturbations which will excite the atom or ion 
into other states.) The novel aspect of a quantum system is that where a classical bit must either 
be a or a 1, a quantum bit can be any suerposition 

|^ = C |0) + C 1 |1), (1) 

where Co and C\ are arbitrary complex numbers, subject to a normalization condition such as 
| Co | 2 + | Ci | 2 = 1, expressing the fact that the system has unit total probability for being in some 
state or other. (We shall leave this normalization out below. This keeps the formulas looking neat, 
and it is easy to reimpose at any point in the argument.) 

The state (|j) has no classical interpretation. In particular it is simply incorrect to say that the 
system is in one of the states |0) or |1) with certain probabilities, in the way that we can say that a 
coin which we toss with our eyes closed is either heads or tails even though we don't know exactly 
which one it is. Such an interpretation of the state leads to wrong predictions for the outcome 
of many measurements. Similarly, saying that the system is simultaneously in states |0) and |1) 
has only limited usefulness. One simply has to accept |TJ) as the state of the system, and use the 
quantum mechanical assignment of measurables to operators to make predictions. 

Next, let us consider a quantum register composed of two two-state systems. The states of this 
register can be written as 0), 1), |1 0) and |1 1). Here |1 0) means that the first two state system 
is in state 1 and the second in state 0. These four states can be regarded as representations of the 
four binary numbers 00, 01, 10, and 11. Again, as in Eq. (|l|), our two two-state-system register can 
be in a superposition 

|^) = Coo|00) + Coi|01) + C 10 |10) + Ch|11). (2) 

Generalizing to a register of N two-state systems, it is obvious that one can represent 2^ numbers 
in this way. The general superposition can now be written as 

2 N -1 

|V) = ]T C n \n), (3) 

n=0 



where we use the notation |n) to represent a state in which the binary representation of the number 
n yields the individual state of each two-state system. 

Suppose now that we wish to evaluate an integer valued function f(n) of the integers. (We 
suppose that both n and f(n) are bounded so that we can represent them with registers of managable 
size.) Let us assume that we have found a means of dynamically evolving our system in such a way 
that the corresponding unitary operator Uf transforms the state \n) into \f(n)}. More precisely, it 
is necessary to consider a computer with two registers, which we call the input and the output. The 
initial state of our computer is taken to be 

^initial = KJO) ou f (4) 

Then the unitary operator Uf is supposed to be such that 

l^) final = ^/l^)i n iti a I=k)inl/W>ouf (5) 

Note that the input register is left unchanged, while the output register contains the value of the 
function. Whether or not the operator Uf can be constructed depends on the function / under 
consideration, and is a separate question, belonging to the theory of algorithms and computability. 

The core concept in quantum computing is to take advantage of the superposition principle. It 
follows from the linearity of the unitary time evolution operator Uf that 

^/E C »inl°>out =E C »inl/(")>ouf (6) 

n n 

Equation ([3J) contains the essence of quantum parallelism. Since the effort required to create or 
implement the evolution Uf is separate and independent of the state of the quantum computer, it 
follows that if we can prepare the initial state of the quantum computer in the superposition on the 
left hand side of Eq. (||), the final state will contain, albeit as a superposition, the function /(n) for 
all values of n. We have performed, at first sight, 2 N function evaluations at one shot. 



1.2 The Read-Out Problem and Shor's Breakthrough 

The parallelism of Eq. (^) is somewhat hollow, however. To obtain information on the function 
values f(n), one must measure the state of the input and output registers. Such a measurement 
can only yield one value of f(n). We have no way of knowing the other function values. Worse, 
the value of n for which the function is evaluated is itself random: the probability of getting n is 
\C n \ . Thus one could not even generate a table of f(n) by repeated runs of our machine without 
substantial redundancy. 

The above difficulty, also known as the read-out problem, stymied the subject of quantum com- 
puters for the decade after its genesis. Interest in the field was primarily in terms of the theory of 
computability and complexity,^, and the thermodynamics of computation, specifically the issues 
of reversibility and energy consumption.]^, ^ The breakthrough came in 1994 with the discovery by 
Coppersmith and Shor of a quantum algorithm for factorization of a composite number with two 
prime factors. jlO| [ll| The key realization is that for certain problems what one needs is not the 
individual function values f(n), but some global property of the function. In the factorization prob- 
lem, the property which is exploited is the period of a periodic function. This can be found using 
Coppersmith's quantum Fourier transform. IK)! We refer readers to a recent review of the number 
theory underlying these algorithms,^) but briefly and heuristically speaking, quantum mechanics 
is intrinsically well set up for the efficient addition of terms multiplied by phase factors required in 
finding a Fourier transform. (Witness the closely related examples of Fraunhofer and iV-slit diffrac- 
tion in optics.) In this way, one can indeed exploit quantum parallelism to do certain calculations 
very effectively on a quantum computer. For a number with L decimal digits, Shor's algorithm 
factorizes it in order L 3 steps.Q The best currently known classical algorithms, by contrast, require 

1 For completeness, it should be noted that Shor's algorithm is probabilistic, in that it requires choosing a random 
number for its operation. The number theoretic method employed by the algorithm fails to give a useful result for 
certain choices of this random number. However, it can be shown that the probability of success can be made to 
approach arbitrarily close to unity with O(L) choices of the random number. Thus this consideration does not alter 
the polynomial time nature of Shor's method. 



Table 1: Truth Table for the Controlled-Not Gate 



Input 


Output 


|0 0) 


|0 0) 


|01) 


|01) 


|10) 


111) 


111) 


|10) 



~ exp[cL 1 / 3 (ln L) 2//3 ] steps, where c is a constant of order unity. Needless to say, this discovery has 
galvanized the subject. The difficulty of factorization is the basis for much modern day cryptogra- 
phy, and the elimination of this difficulty would have obvious repurcussions. In addition, it has been 
surmised that a generic quantum computer may be used for the simulation of a general many body 
quantum mechanical problem. || 

1.3 Primitive Quantum Gates 

It remains to say something about how the unitary operator Uf is to be constructed. One would 
like to do so by composing a small number of basic or primitive operations operating on a small 
number of bits at a time. Indeed, a one-bit gate consisting of a rotation, and a two-bit gate called 
the controlled-not (or C-NOT) suffice to build any operator Uf, in the same way that an AND and 
a NOT suffice for a classical computer. (Equivalently, one could say that any unitary matrix acting 
on an n dimensional space can be built by multiplying n(n — l)/2 unitary matrices that act on two- 
dimensional subspaces at a time. The general point that such simple operations can be composed 
to form arbitrarily complex operators is obvious to any one who has ever solved linear equations 
by Gaussian elimination or solved matrix problems numerically, but readers who wish to see this 
codified as theorems should consult Refs. |(| [b| [l5| . The physical point which is important is 
that it is necessary to have a means of introducing correlations amongst the states of two different 
two-state systems, and this is done via the controlled-not. Any other two-bit gate that was not too 
close to the identity operation would work just as well.) The one-bit rotate acts on a single bit as 
follows: 

£70X0)10) +Ci|l» = (Cocos^ + Cisin^lO) + (-C sin</> + C x cos$|l), (7) 

while the C-NOT acts on any pair of bits as per the truth table in Table 1. The action of this gate 
is best described in words. It leaves the first (or control) bit unchanged, while the second (or target) 
bit is flipped from 1 to or vice versa if the control bit is a 1 , and left unchanged if the control bit is 
a 0. It should of course be remembered that this operation also acts on linear superpositions. Thus, 

o|01) + 6|ll) C ~ NOT > o|01) + 6|10>. (8) 

It should be noted in this connection that the C-NOT is defined with respect to a very definite 
choice of basis states, which is referred to as the computational basis. It should also be noted that in 
practical implementations, it may be more useful to construct a few more gates than just the C-NOT 
and the one-bit rotate. Examples of such gates and simple circuits like adders and multipliers may 
be found in Refs. § g. 

The above discussion of the construction of the evolution operator Uf is very general, and does 
not say how the individual gate operations are to be carried out. There are at least two broadly 
different types of quantum computers that have been discussed in the literature. In Feynman's 
original conception^, the idea is akin to having an array of sites between which an electron can 
move, and to arrange the terms in the Hamiltonian to be such that matrix elements between different 
sites implement the desired gates. It is important that the connections be local in order that the 
gates remain simple. A computation then proceeds in analogy with an electron wave packet coursing 
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Figure 1: Schematic of the ion-trap quantum computer. 



through the sites. The second conception, which is the one more commonly discussed today is to 
have a temporal sequence of gate operations applied through external perturbations to a fixed set 
of quantum bits. 

With this somewhat long preamble, we now turn to the physical implementation of a quantum 
computer, especially that based on ion traps. 

2 THE ION-TRAP QUANTUM COMPUTER 

The desiderata in making a real quantum computer are almost too obvious to state. One needs 
first and foremost to have a means of controlling and driving the bits individually. Secondly, one 
needs to have low dissipation, or interaction with the environment. Equation (|^) assumes that the 
computer forms a closed system, isolated from the rest of the world. In reality, interactions with 
the environment will end up entangling the state of the computer with that of the environment. 
These interactions are most likely to be such that the correlations amongst the computational 
degrees of freedom (or quantum bits) get lost over time. This point is well known from the study of 
problems like NMR, where also one has many quantum systems, each with a Hilbert space of small 
dimensionality, but interacting with an environment — the lattice, motional degrees of freedom of 
the atoms, electromagnetic radiation, etc. 

Even after the above conditions are met, it is not easy to implement a good two-bit gate. In 
addition to the ion-trap proposal^] that we shall describe in more detail, another one based on caity 
QED[fl6| seems promising. In fact, experimental two-bit gates, inspired by each of these schemes, 
were demonstrated shortly after their proposal. |l7], |l8) 

The ion-trap QC utilizes an array of N identical ions in a linear rf Paul trap, each of which 
can be independently addressed and driven by a laser. (See Fig. 1.) The ions are cooled to nearly 
zero temperature. Such arrays have been studied for some time now because of the opportunities 
they offer for spectroscopy, and developing frequency standards. [jl9|, |2(i[| The ions have a nonuniform 
spacing determined by the Coulomb repulsion between the ions and the effective trap potential. A 
quantum bit is composed of two internal states of an ion, which we denote by |e) and \g), and one-bit 
gates are executed by applying suitable pulses (ir/2, it, etc.) to the ions. The clever idea which 
enables the execution of a two-bit gate is to use the center of mass motion of the entire array as 
an additional computational degree of freedom. Let wo be the g <-> e transition frequency of the 
ions, and let lo z be the center of mass vibrational frequency. A laser pulse with carrier frequency 
ujq ± lj z acting on the ith ion entangles the internal state of this ion with that of the center of mass 
mode. As discussed in Ref. ||, a sequence of three such pulses, the first and third applied to ion i, 
and the second to ion j, have the net effect of a two-bit gate in which ions i and j are the control 
and target bit, respectively. In fact, by also applying one-bit rotate operations to ion j, this gate 
can be exactly turned into the C-NOT. A sequence of one- and two-bit gate operations can then be 
executed, yielding any desired computation, as discussed in Sec. [j]. 



2.1 Spontaneous Emission Decoherence 



The advantages of the ion-trap QC are its conceptual simplicity, and high degree of isolation from 
the environment. In fact, the proposal is regarded seriously enough that an effort is underway to 
build a proof-of-principle prototype. |2l| The most obvious source of decoherence is the spontaneous 
emission from the ions. It seems reasonable that effects like superradiance, in which the proximity 
of the ions to one another affects this decay, can be neglected, and that each ion can be treated 
separately. If we assume that half the ions are in an excited state at any given time on the average, 
then in a computer with N ions, any computation must be completed in a time of order r ra d = 2t s /N, 
since even one spontaneous emission event destroys the phase coherence among the computational 
basis states of the QC as a whole. One can minimize this source of decoherence by employing ions 
where the |e) state can not decay to any state other than \g) and the decay in question is dipole 
forbidden. For example, in Ba + ions, the decay between the 5d 2 D 5 / 2 multiplet (one member of 
which serves as |e)) and the 6s 2 S , 1 / 12 multiplet (which supplies the \g) state) is via an E2 process, 
for which r s = 30-70 s.[] The drawback, as noted by Plenio and Knight[^2| is, that by the general 
relation between Einstein A and B coefficients, the coupling of the laser to the g «-> e transition is 
also weak, and the time required for all gates is correspondingly long. One can not shorten this time 
without limit by increasing the laser power without starting to excite the ion into higher levels, or 
even ionize it further, via two-photon absorption. The radiative decoherence can be decreased by 
working with a A system and using Raman pulses, but it can not be totally eliminated. |2^| 



2.2 Vibrational Decoherence 

The second source of decoherence in the ion-trap QC is the ionic vibration. To understand this, 
suppose (see Fig. 1), ion j is displaced from its equilibrium position. This creates an excess electric 
field or electric field gradient at any other ion. say i. The time evolution in the |e), \g) space of 
ion i is thereby altered, and in the long run, the QC is not in the intended state, i.e., we have 
decoherence. Let us denote the decoherence time due to this process by T v ;b. A proper calculation 
of this time is not so easy as that of r ra d as all the ions are coupled to the same bath, and can not 
be treated individually. A naive consideration based on only a few ions leads to an overestimate, 
and is misleading. The calculation presented in the Appendix shows that scales more rapidly 
than A, where N is the number of ions. More precisely, we find, 

t~} ~ A 1/2 , (9) 

where Q is a quadrupole transition matrix element (appropriate to the Ba + example), ojt is a typical 
transverse normal mode frequency for the ion array, q and m are the ionic charge and mass, and so 
is the minimum spacing between the ions which occurs at the center of the array. We can write t~^£ 
more explicitly as a rate by noting that q 2 — tow 2 dp, where do is the the trap length scale parameter 
[see Eq. (|3|)], and that 

Q 2 oc ft/r.fcg, (10) 

where k = oj / c. It follows that 

1 N 1 ' 2 /d \ 3 uj 2 1 



Tvib T s \S J LU LO t (fc So) 



(11) 



2.3 Discussion of Results 

Equation ( pd| ) has interesting implications for the scaling of r V ib with N. We can not naively take this 
as A 1 / 2 because the behaviour of u x , uJt, and so as A is increased, depends on how the trap operating 
conditions are varied. Suppose so is held fixed as A is increased. Then (rfo/so) 3 ~ N 2 /\nN and 
T~ ib ~ A 5 / 2 / In A. In this case, however, the longitudinal voltage on the trap electrodes, which is 

2 Other ions in groups IIA and IIB, such as Ca+ and Hg+, have similar level schemes, and long spontaneous decay 
times. 



proportional to w 2 , varies as (\nN)/N 2 . Since the time needed for a two-bit gate varies as uj^ 1 , 
the total computational time goes up. Secondly, the longitudinal confinement becomes weaker, 
and non-linearities in the trapping potential, and electrode patch voltages become more important. 
Suppose, on the other hand, that the trap voltages, and therefore, uj z and u>t, are held fixed as N 
increases. Then tT^ ~ A 35 / 6 (ln A)~ 8 / 3 . [The variation is as A 9 / 2 (ln7V)~ 2 for an El transition.] 
Now, however, the minimum inter-ion spacing Sq varies as ~ A~ 2 / 3 . This may make it difficult to 
optically resolve and address individual ions, which is basic to the operation of the QC. It seems 
likely that some compromise between these two extremes will have to be sought, depending on 
engineering considerations. A general point is worth noting in this connection. Since the radiative 
and vibrational decoherence processes are independent, the total decoherence time of the computer 
is given by adding their rates: 



td is the useful window of time in which any computation must be finished. Thus if one of the two 
rates turns out to be much larger than the other, we can relax the design considerations on the 
smaller rate, and focus on ways to reduce the larger one. 

Let us estimate T v ;b using the example of Ba + ions and the levels mentioned above. The frequency 
uj a = (2vr)1.7 x 10 14 Hz. We take uj z /2tt = 100 kHz, and w t /27r = 20 MHz. This yields d Q = Ufxm. 
For N — 1000, we obtain r v ;b — 10 4 r s , which is surprisingly large. (It is even larger in comparison 
to T ra d = t s /N .) The drawback is that Sq ~ 0.5 /im with the same parameters. This runs into the 
difficulty with optical resolution mentioned above. We have not explored compromise variations of 
trap parameters in detail. 

It is natural to ask if we should not have anticipated r ra d being so much larger than r V ib- The 
non-trivial scaling with N and the fact that the relationship would reverse for larger N makes us 
believe that the answer is no. 

2.4 Conclusion 

In summary, vibrational decoherence is not a significant problem in the ion-trap QC for N < 10 3 , 
as originally envisaged by Cirac and Zoller.j^] This optimistic result should be tempered somewhat 
however. The technical difficulties in working with 10 3 ions are enormous, and the spontaneous 
emission decoherence itself is nothing to sneeze at. It is possible that the A system with Raman 
transitions will mitigate this difficulty, but this needs to be explored more fully. A qualitative 
argument suggests that the adiabatic suppression of t~^£ will not be as effective in this case, however. 
If this turns out to be true, it would provide a nice example of the compromise between different 
types of decoherence discussed above. 



A CALCULATION OF VIBRATIONAL DECOHERENCE 



We give here an approximate calculation of the vibrational decay time. The calculation is based on 
Ref. (J] which should be consulted for a more precise description. 

We denote the equilibrium axial position of the j th ion by Zj , and its deviation from equilbrium 
by Uj . We will consider ijne transition of electric quadrupole (E2) type, and denote the ionic 
charge by q, and the quadrupole transition matrix element by Q, ignoring all vector and tensor 
indices, here and in what follows. The coupling between the vibrations and the eg space of the ith 
ion is described the Hamiltonian 




(12) 



TIME 




(13) 



The effect of this perturbation is small for two reasons. First, its magnitude is small, \Vi\ -C %uiq. 
Second, the time scale for its variation is determined by the ion array's normal mode frequencies, 



which are of order 100 MHz at most, and thus much less than luq, which is an optical transition 
frequency of order 100 THz or more. Formally, 



Yi 

Vi 



< wo- 



(14) 



The perturbation is thus slow, and the internal state of the ions can follow the vibrational distur- 
bances adiabatically. If this following were perfect, we would have no real excitations or transfer of 
energy. Decoherence arises solely because of corrections to the adiabatic approximation, and we can 
thus expect the final answer for r^ b to contain a small factor like (u> v /cJo) a where a is a positive 
exponent, and u v is a typical vibrational time. We now turn to estimating this effect. 

A.l Adiabatic Approximation 

Let us map each two-state ion onto an equivalent spin-1/2 system, with |e) and \g) being the up 
and down spin states. The ith spin sees effective magnetic fields B z = Hloq and B±(t) = Vi(i) (see 
Fig. 2), where we take Vi{t) to be a specified time dependent c-number perturbation. The total 
Hamiltonian for this ion can thus be written as 



Hi = ^huj <J z + hf(t) ■ a, 



(15) 



where f(t) = (B x ,B y> 0)/2h, and the tr's are the usual Pauli matrices in the equivalent spin space. 
Note that f has no z component. 

We wish to study the time evolution of the spin for a general f(t) which is small and slow. To 
this end, let us write a general spin state as 

\i/>(t)) = u + (t)e- luJ0t/2 \+) + U _(i)e^ 0t/2 |->- (16) 

where |±) denotes the eigenstates of S z with eigenvalues ±1/2. Schrodinger's equation gives 

iu± = e ±i »° t Mt)u T (t), (17) 

where /± = / x ± ify Since f± varies very slowly, we try a solution to Eq. ( p~T| ) of the form 

u±(t) =a±(t)+p±(t), (18) 

where a± and [3± are functions that vary rapidly and slowly, respectively. The latter is almost a 
constant over a period 2tt/u)o, and the former almost averages to zero. 
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Figure 2: Equivalent magnetic field description of ith ion. 



We now substitute Eq. ( |l8| ) in Eq. ([l?]), and separately equate the fast and slowly varying parts. 
This yields 



zd± = e^'/T^/MO, (19) 
i$± = e ± ^ t / T (t)« T (t). (20) 

The next step is to integrate these equations. Consider Eq. ( |l9| ) first. To a first approximation, we 
can regard f± and (3± as constants. Integration then gives 

a± = T^ " 1 e ±ia;ot / T (^W- (21) 
Substitution of this result in Eq. ( p(i| ) yields 

\{(t)\ 2 

z/3± = ±JiWL^(t). (22) 

Integrating this equation, we obtain 

/3±(t)=exp(T**(t))jft t (0), (23) 

where t 

$( t ) = / tft' M^IL. (24) 
Jo 

We next note that since by Eq. (pl|), |a±//3±| ~ |/±/ w o| Cl,«±~ /3± in leading order in f±/cJo. 
Equations @, @, and @ thus give us the desired expression for the time dependence of a 
general state. 

The above equations also give us the amount of decoherence. Suppose the initial state of the 
spin is 2 _1 / 2 (|+) + |— )), i.e., u±(0) = 2~ 1 / 2 . (A product of states of this type for each quantum bit 
is the starting point in several quantum algorithms.) Let the state that would be obtained at time 
t in the absence of the perturbation f be \ipo(t)). The states \ipo(t)) and \i/>(t)) are the states of the 
quantum bit in the ideal and actual QC, without and with decoherence. The decoherence is given 
by the overlap 

(Vo(t)|V(*)>=cos($(t)). (25) 

Another way to understand these results is that because \B±\ -C B z , the precession axis for the 
spin can be taken to be z at all times to very good approximation, and the instantaneous precession 
frequency can be taken as 

<4 = H + ^ 2 A 2 ) 1/2 « + (26) 

In ujq 

Hence, the time dependence of the states |±) is given by exp (±i J* dt' u>' 0i (t')/2\ |±). [Note that 
$(t) is just half the difference between J Q u)' oi (t') dt' and to^t.] 



A. 2 Single and Many Ion Decoherence Times 

The decoherence time for the ith ion is approximately equal to the time required for its phase to 
deviate from its unperturbed value by 7r: 7~i ~ ir/(ui 0i — lu ). Using Eqs. (|l^) and (p6|), we obtain 

where the angular brackets denote some kind of average. In principle this average should be calcu- 
lated by expressing the displacements u in terms of the normal modes of the array, which should be 
described by a thermal density matrix at the temperature T to which they have been cooled. (Ad- 
ditional correlations from the center of mass mode excitations inherent in the two-bit gates should 
be excluded since they are part of the intended time evolution and do not generate decoherence. 



We have already left this part out in deriving Eq. ( p7| ) by taking f to be slow.) We are interested 
in obtaining an order of magnitude answer. Therefore, we will ignore the normal mode structure, 
and assume each ion to vibrate independently with some average frequency. Secondly, to obtain the 
best-case answer, we will assume that T <C hu> z /kB, and simply set T — 0. To be more concrete, we 
write 

(ujU k ) = 5jk, (28) 

moj t 

where u)t is a typical transverse mode frequency, and m is the ionic mass. It follows that 

, q 2 Q 2 ^ 1 



lixhmLOQUit ^-^ (zi — Zj) 

There are two reasons for using a transverse mode frequency in Eq. (|2^). First, in the linear Paul 
trap, these frequencies are generally higher than the longitudinal ones. Second, one can arrange via 
a J z selection rule for the longitudinal modes not to excite any j <-> e transitions. 

We now combine the individual ion decoherence times to obtain r v ib for the QC as a whole. The 
obvious procedure of adding the rates is not quite correct. The reason can be seen from Eq. (p5|). The 
overlap between the actual and intended states of the zth spin is more like cos(i/rj) than exp(— t/ri). 
Multiplying these overlaps for all spins, we obtain P(t) ~ Hi cos 2 (t / Ti) ~ exp(— t 2 /r 2 ih ) with 



E-r 2 - 00) 



A more careful justification for this result can be found in Ref. 0. We note here that if we simply 
add t" 1 to obtain t^c, we overestimate tT 1 by a multiplicative factor of TV 1 / 2 ; if this were done, 



the N 1/2 in Eq. (ITT]) would change to N. 



A. 3 Continuum Approximation for Ion Array 

Equations ( p9| ) and ( |30| ) provide us with a formal answer for the vibrational decoherence rate. To 
obtain a more useful result, however, it is necessary to perform the sums over the lattice positions. 
For small N, say 25 or less, these sums are best done numerically, but for larger N, we can evaluate 
them by approximating the ion array as a continuum. 

Let us define s(zi) to be the average of the local spacing between the ion at Zi and its two nearest 
neighbours. We expect that for large N, this spacing will vary slowly as we move along the array. 
Hence, it is a good approximation to treat z as a continuous variable, and on this basis we will find 
an expression for s(z). We will also find the total length of the array as a function of N. 

Consider an ion at position z, and let us denote the separations to its immediate neighbours to 
the left and right by s_ and s + respectively. The Coulomb force on the ion at z from these near 
neighbours is given by 

i^n° Ul = 9 2 (s= 2 - s+ 2 ) « 2q 2 S - 2 (ds/dz), (31) 

since s+ — s_ « s(z)(ds/dz). For ions not too close to the edge of the array, the distances to the 
second, third, fourth, neighbour pairs are approximately doubled, tripled, quadrupled and so on. 
The forces from these successively distant neighbour pairs are smaller than i 7 ^ " 1 by factors of 4, 9, 
16, etc. Since, 53n=i n ~ 2 = ^V^' and since the sum converges very rapidly, we can write the net 
Coulomb force on the ion at z as 

^ct ul = (7T 2 q 2 /3 S 2 )(d S /dz). (32) 

Equating this to the opposing spring force muj 2 z from the trapping potential, we obtain a differential 
equation for the spacing function: 

7T 2 ds Z 

2^{z)Tz = ■ (33) 

We have introduced 

d Q = (q'/mcol) 1 / 3 (34) 



as a natural length scale for the trap. (One can easily show that the ion spacing is of order do for 2 
or 3 ions in the trap.) 

To integrate Eq. (|33|), we define the total length of the array to be 2L. Choosing the center to 
be at z = 0, we obtain 

<L) ~W) = ~2^4 {L2 ~ z2) - (35) 

To make use of this result, however, we need expressions for L and s(L). To find s(L), we use same 
argument as was used to obtain Eq. (|32"|), although the approximation is clearly not as good now. 
We assume that the spacing between the ions near the end of the chain is uniform and equal to s(L). 
The Coulomb force on the last ion is then ir 2 q 2 /6s 2 (L). Balancing this with the spring force muj 2 , 
we obtain s(L) w n(d 3 ) /6L) 1 ^ 2 . Even if this argument is not watertight, it shows that we can ignore 
s~ 1 (L) compared to 3L 2 /2ir 2 d,Q in Eq. (p5|). This yields 

S (z) = s (l-z 2 /L 2 )-\ (36) 

where we have introduced the minimum ion spacing (attained at z = 0): 

s = s (0) = 2TT 2 dl/L 2 . (37) 

It still remains to find L(N). If we denote the ion number at position z by n(z\, then dn/dz — 
\/s(z) in the continuum approximation. Integration of this result along with Eq. (|36| ) gives 

L = d^ 2 N/2f/\ (38) 
s = 6.81d iV~ 2/3 . (39) 

One can also show that the mean spacing varies as N~ 2 ^ 3 In AT. 

The treatment of the ends of the chain above is not fully satisfactory. Another continuum 
approach is due to Dubin.[^3| He regards the array as fluid of total charge qN. It is known that in a 
harmonic potential such as that of the trap, such a fluid forms an ellipsoid of revolution of uniform 
charge density. When the trap is much stiffcr in the transverse than the longitudinal direction, i.e., 
u>t 3> u) z , the ellipsoid has a total volume AnNd^. Equating the semi major axis of the ellipsoid to 
L, we find the semi minor axis to be i^Nd^/ L) 1 / 2 . The spacing s{z) is now given by geometry. Let 
A(z) be the cross sectional area of the ellipsoid at an axial distance z from the center. The volume 
s(z)A(z) clearly contains one unit q of charge, from which it follows that 

This agrees with Eq. ((3^), but so has a different form. To complete the solution, we must find 
this form, or equivalently, L(N). The fluid model answer for L depends on to z /uJt, which is clearly 
wrong if the linear structure is stable. Dubin therefore resorts to a local density functional theory 
to estimate the discreteness correction to the Coulomb energy. He then minimizes the sum of this 
correction, the fluid drop self energy, and the trapping potential energy. The resulting array length 
is independent of uj z /idt and is given by 

L 3 = 3Nln(c N)dl, (41) 

where Co = 6e 7 ~ 13 / 5 « 0.8, and 7 is Euler's constant. We also obtain 

s = 4L/3N = 1.92A^ 2/3 [ln(0.8A0] 1/3 d o . (42) 



These results differ from Eqs. (38) and (|39j) logarithmically in N. 

It is interesting to compare these results with some recent semi-numerical work by Meyrath and 
James. p4| They integrate Eq. (|36"|), write the result as n(z) = az — bz 3 , and fit a and b to power laws 
in A^ instead of trying to relate them to so and L. They then invert this cubic equation for z(n), 
and find that the answer agrees quite well with numerics for A^ > 25. Likewise, Eq. ( |42] ) provides a 
very good fit to the numerical results. 



A. 4 Lattice Sums 



Equations (|40|442|) can be used to find the sums in Eqs. (29) and (BCt). The first type of sum, 



can be very simply evaluated as 



^z)«2 S »(zO£4 = #T ( 44 ) 

S n (Zi) 

This result should hold well for all i except very close to the ends, since the exponent n is big (3 for 
an El decay, 4 for E2). 

The second type of sum is T n = s~ n (zi). Writing Ai w dz/s(z), we can approximate it by 
an integral: 



dz 

-L 

Substituting Eq. ( |36| ) and performing the integration, we obtain 

1/2 



Tn^-^ll-r^) , (46) 



s™ +1 V 4n + 7 



where we have also used an asymptotic formula for (3(n + 2, 1/2). 

The combination of Eqs. @ and © with Eqs. © and © leads to the results ©-([□]) 
quoted in Sec. 2. 
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